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Interacting electrons in a semiconductor quantum dot at strong magnetic fields exhibit a rich set 
of states, including correlated quantum fluids and crystallites of various symmetries. We develop in 
this paper a perturbative scheme based on the correlated basis functions of the composite-fermion 
theory, that allows a systematic improvement of the wave functions and the energies for low-lying 
eigenstates. For a test of the method, we study systems for which exact results are known, and 
find that practically exact answers are obtained for the ground state wave function, ground state 
energy, excitation gap, and the pair correlation function. We show how the perturbative scheme 
helps resolve the subtle physics of competing orders in certain anomalous cases. 
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There is a strong motivation for developing theoretical 
tools for obtaining a precise quantitative description of 
interacting electrons in confined geometries, for example 
in a semiconductor quantum dot, because of their possi- 
ble relevance to future technology |l| . Exact diagonaliza- 
tion is possible in some limits but restricted to very small 
numbers of electrons, and does not give insight into the 
underlying physics. For larger systems, one must nec- 
essarily appeal to approximate methods. The standard 
Hartree-Fock or density functional type methods provide 
useful insight, but are often not very accurate for these 
strongly correlated systems. The aim of this paper is 
to demonstrate that a practically exact quantitative de- 
scription is possible for a model quantum dot system, 
facilitated by the ability to construct low-energy corre- 
lated basis functions. 

Our concern will be with the solution of 
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which contains N interacting electrons in two dimensions, 
confined by a parabolic potential and subjected to a mag- 
netic field. The parameter rat, is the band mass of the 
electron, c^o is a measure of the strength of the confine- 
ment, e is the dielectric constant of the host semiconduc- 
tor, and rjk = \rj — r^|. We will consider the limit of a 
large magnetic field (u) c = eB/mi,c 3> ujq), when it is a 
good approximation to take electrons to be confined to 
the lowest Landau level (LL). In that limit, the energy 
eigenvalues have the form E(L) = E C (L) + V(L) where 
the contribution from the confinement potential is explic- 
itly known as a function of the total angular momentum 
L: E C (L) = (h/2)[Q - lu c ]L, with fi 2 = u? c + 4lu%, and 
V(L) is the interaction energy of electrons without con- 
finement, but with the magnetic length replaced by an 
effective magnetic length given by I = ^Jh/m^Vl. Thus, 
the problem is reduced to finding the interaction energy 
V (which will be quoted below in units of e 2 /e£) as a func- 
tion of the angular momentum L. Exact results, known 



for a range of N and L values from a numerical diag- 
onalization of the Hamiltonian, provide a rigorous and 
unbiased benchmark for any theoretical approach. Exact 
studies have shown [2| a correlated liquid like state at 
small L, but a crystallite at relatively large L, as may 
be expected from the fact that the system becomes more 
and more classical as L increases. (The ground state in 
the classical limit is a crystal |3|.) Hartree-Fock studies 
have been performed for the quantum crystallite |J, |5(. 

We will apply the composite fermion (CF) theory [6j 
to quantum dot states 0, S @ ■ The central idea is a 
mapping between strongly interacting electrons at an- 
gular momentum L and weakly interacting electrons at 
L* = L — pN(N — 1). In particular, a correlated basis 
for the low energy states of interacting electrons 
at L can be constructed from the trivial, orthonormal 
Slater determinant basis for non-interacting electrons at 
L* , denoted by }, in the following manner: 
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oj — .xj j — iyj denotes the position of the jth elec- 



Here, zj = Xj 

tron, 2p is the vorticity of composite fermions, and V 
indicates projection into the lowest LL. (Electrons at L* 
in general occupy several Landau levels.) The symbol 
a = 1,2, ••♦,£)* labels the D* Slater determinants in- 
cluded in the study. In general, the basis {^/^} is not 
linearly independent, so its dimension, -Dcf, may not be 
equal to D* (Dcf < D*). The advantage of working 
with the correlated CF basis is that £>cf is drastically 
smaller than D ex , the dimension of the lowest LL Fock 
space at L (which is also the dimension of the matrix 
that must be diagonalized for obtaining exact results). 
Fig. illustrates some basis functions at L = 95. 

At a back-of-the-envelope level, one can compare the 
exact interaction energy at L to the kinetic energy of free 
fermions at L* , with the cyclotron energy treated as an 
adjustable parameter 0- That reproduces the qualita- 
tive behavior for the L dependence of the exact energy 
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FIG. 1: Schematic depiction of Slater determinant basis states 
for N = 6 electrons at L = 95, which maps into L* = 5 with 
2p = 6. The single electron orbitals at L* = 5 are labeled by 
two quantum numbers, the LL index n — 0, 1, and the an- 
gular momentum / = — n, — n+1, ... The x-axis labels n + l and 
the y-axis n. The dots show the occupied orbitals forming the 
Slater determinants =5 relevant up to the first order. The 
state shown at the top left has the lowest kinetic energy (if the 
kinetic energy is measured relative to the lowest Landau level, 
then, in units of the cyclotron energy, the total kinetic energy 
of this state is two). The other nine states have one higher 
unit of kinetic energy. The basis states ^>^ =95 are obtained ac- 
cording to Eq. through multiplication by Ylj<k( z J ~ z kY\ 
which converts electrons into composite fermions carrying six 
vortices. That is shown schematically by six arrows on each 
dot. The single state at the top is relevant at the zeroth order, 
and all ten basis states are employed at the first order. (In 
fact, there are a total of 12 linearly independent states <3>^ a t 
the first order for L* — 5, but they produce only ten linearly 
independent states at L = 95.) 



for small L [7J , but discrepancies are known to appear at 
larger L H0]. 

For a more substantive test of the theory, it is neces- 
sary to obtain the energy spectrum by diagonalizing the 
Hamiltonian of Eq. 1 in the correlated basis functions of 
Eq. J2J. The CF-quasi-Landau level mixing is treated 
as a small parameter, and completely suppressed at the 
simplest approximation, which we refer to as the zeroth 
order approximation. Here, the correlated basis states at 
L are obtained by restricting the basis {$ Q } to all states 
with the lowest kinetic energy at L* (with p always cho- 
sen so as to give the smallest dimension). 

Diagonalization in the correlated CF basis is techni- 
cally involved, but efficient methods for generating the 
basis functions as well as all of the matrix elements re- 
quired for Gram-Schmidt orthogonalization and diago- 
nalization have been developed using Metropolis Monte 
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TABLE I: Exact ground state energy (14 x ) and the ground 
state energy obtained from the zeroth (V^p) and the first- 
order (Vq]) ) CF theory for N — 6. The dimensions of the 
bases diagonalized are D cx , DjSp and -Dq F , respectively. The 
statistical uncertainty arising from Monte Carlo sampling is 
given in parentheses. 



Carlo sampling. We refer the reader to earlier literature 
for full details [I HfH . 

A diagonalization of the Hamiltonian in the zeroth 
level basis produces energies and wave functions for D^ F 
low-lying states. The interaction energy and the wave 
function for the ground state will be denoted V^p and 

^q F , respectively. We have carried out |Tl| extensive 
calculations for a large range of L for up to ten particles, 
and found that the CF theory reproduces the qualita- 
tive behavior of the energy as a function of L all the 
way to the largest L for which exact results are known. 
We show in Table [i] results for N = 6 electrons in the 
angular momentum range 79 < L < 108, which spans 
both liquid and crystal-like ground states. Ref. 01 
shows a comparison of the exact pair correlation func- 
tion with that calculated from 

^cf for L = 95 (where 
there is a unique CF wave function); surprisingly, the 
CF theory, originally intended for the liquid state, au- 
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TABLE II: Overlaps between exact ground states and CF 
ground states obtained at the zeroth (O' - 1 ) and the first 
order(O^). The statistical uncertainty from Monte Carlo 
sampling does not affect the first three significant figures. 
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FIG. 2: Pair correlation function for N — 6 electrons at 
L = 99. The position of one particle is fixed on the outer ring, 
coincident with the position of the missing peak. The ground 
state wave function used in the calculation is obtained from 
(a) exact diagonalization; (b) the zeroth-order CF theory; (c) 
the first-order CF theory; (d) the second-order CF theory. 
The "noise" in (a) and (d) results from the relatively large 
statistical uncertainty in Monte Carlo because of the more 
complicated wave function. 



tomatically produces also a crystallite at large L, even 
though no crystal structure has been put into the the- 
ory by hand [lj- Table HD gives the overlaps defined as: 
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While the zeroth level description is quite good, the fol- 
lowing deviations from the exact solution may be noted, 
(i) The overlaps are in the range 0.70-0.94, which are 
not as high as the overlaps (~0.99) for incompressible 
ground states in the spherical geometry, (ii) The ener- 
gies are within 0.5% of the exact ones, which is quite 
good but could be further improved, (iii) In the crys- 
tallite, the particles are somewhat less strongly localized 
in the CF wave function than in the exact ground state. 
See Ref. [ll| . (iv) The symmetry of the crystallite is pre- 



dicted correctly with the exception of L = 99, where the 
CF theory predicts a (0,6) crystallite [Fig.^b)], that is, 
with all six particles on an outer ring, whereas the exact 
solution shows a (1,5) crystallite [Fig. Eta)], which has 
five particles on the outer ring and one at the center, (v) 
A successful theory must explain not only the ground 
state but also excited states, especially the low-energy 
ones. We have considered the gap between the two low- 
est eigenstates. The zeroth-order theory does not give, 
overall, a satisfactory account of it. In some instances 
(e.g., L = 81,85,90,95,99,105,106 for N = 6), the CF 
theory gives no information on the gap, because the ba- 
sis contains only a single state here (Dq F = 1); in many 
other cases, the gap predicted by the CF theory is off by 
a factor of two to three. 

These discrepancies have motivated us to incorporate 
CF-quasi-LL mixing perturbatively. (We stress that CF- 
quasi-LL mixing implies LL mixing at L*, but the basis 
states at L are, by construction, always within the lowest 
LL.) At the first order, we include basis states <& L at L* 
with one more unit of the kinetic energy, which produces 
a larger basis at L through Eq. The basic idea is 
illustrated for the case of N = 6 and L — 95 (L* = 5, 
2p = 6) in Fig. [2 

In a similar way we have constructed correlated ba- 
sis functions at each angular momentum in the range 
79 < L < 108. As shown in Table EQ-D^, the dimen- 
sion of the basis in the first-order theory is larger than 
Dgp but still far smaller than D cx . A diagonalization of 
the Hamiltonian in this basis produces the ground state 
energy Vq F and ground state wave function Leav- 
ing aside L = 99, which we shall discuss separately, the 
following observations can be made: (i) The energies are 
essentially exact. As shown in TablcQJ the deviation from 
the exact energy is reduced to <0.1%, in fact, to <0.05% 
in most cases, (ii) The overlaps from the first-order 
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given in Table [TT| They are uniformly excellent (0.98- 
0.99) in the entire L range studied, (iii) The improve- 
ment by the first order perturbation theory is also man- 
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FIG. 3: Comparison of the exact excitation gaps (□) for 
N = 6 with the gaps obtained in the first-order CF theory 
(•)• 
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99 86499 76 1.9189 1.9193(3) 0.995 

TABLE III: Comparison of the second order CF theory with 
exact results for the L = 99 ground state. is the dimen- 
sion of the correlated CF basis, Vq F is the CF prediction for 
the ground state energy, and is the overlap between the 
CF and the exact wave functions. 

ifest in the pair-correlation functions, which are now in- 
distinguishable from the exact ones at arbitrary L. That 
is not surprising, given the high overlaps, (iv) As seen 
in Fig.|3]the first-order theory reproduces the qualitative 
behavior of the excitation gap as a function of L, and also 
gives very good quantitative values. The maximum gaps 
are correlated with the states where a downward cusp 
appears in the plot of V(L), consistent with the higher 
stability of these ground states. 

Finally, we discuss the case of L = 99. Here, the ze- 
roth order CF theory predicts a wrong symmetry for the 
crystallite [Fig. Efb)]. As seen in Fig. 0{c), the first or- 
der correction also fails to recover the correct symme- 
try. That is also reflected in the fact that the modified 
ground state of the first-order theory yields a relatively 
small overlap of ~ 0.86, and the energy is off by a rela- 
tively large 0.15%. A closer inspection of the correlations 
in Fig. [2tc) reveals a slight broadening of the hexagonal 
structure in the outer-ring, combined with an appear- 
ance of a small mound at the center, suggesting that 
the structure here may be a superposition of (0,6) and 
(1,5) crystallites. This has motivated us to incorporate 
the next (second) order corrections. The basis dimension 

(2) 

D CF is now further enlarged, but Monte Carlo still pro- 
duces reliable results. The CF ground state wave func- 
tion obtained at this level is extremely accurate: the pair 
correlation function shown in Fig. |2Id) compares well to 
the exact one, and, as seen in Table ITTT1 the energy and 
the overlaps are close to perfect. The origin of the diffi- 
culty can be understood from the fact that the (0,6) and 
the (1^5) crystallites are nearly degenerate in the classical 
limit 3], making them both competitive. 

We have focused in this work on cases where the exact 
results are known, because the aim was to test the appli- 
cability of the CF theory to quantum dots. The theory 
can be extended to much larger systems, where exact di- 
agonalization is not possible; in such cases, one would 
need to increase the accuracy perturbatively until suf- 
ficient convergence is achieved. The method developed 
here should also prove useful for multiple coupled quan- 
tum dots [l3|_and rapidly rotating atomic Bose-Einstein 
condensates 11-41 . 
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